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Abstract 

We have calculated the adsorption energy of carbon monoxide on a mono- 
layer of copper adsorbed on the (111) face of cubic zirconia. We investigate the 
structural parameters of three phases of bulk zirconia (cubic, tetragonal, and 
monoclinic) and find excellent agreement with experiment. We have also an- 
alyzed the structural relaxation of both the stoichiometric and reduced (111) 
surfaces of cubic zirconia (c-Zr02). For adsorption of copper on c-Zr02, we 
find that the preferred binding site is atop the terminal oxygen atom, favored 
by 0.3 eV over other high symmetry sites. We compare CO adsorption on 
zirconia-supported copper to the results of carbon monoxide on copper (100) 
(S. P. Lewis and A. M. Rappe, J. Chem. Phys. 110, 4619, (1999).) and show 



1 



that adsorption on oxide-supported copper is over 0.2 eV more stable than 

adsorption on the bare surface. 
82.65.My, 71.15.Mb, 82.65.Jv 
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I. INTRODUCTION 



Oxide-supported transition metals have received much attention due to their utility in 
automotive catalysis and other catalytic processes. Use of an oxide as a support material 
for transition metals has clear economic benefits because this reduces the amount of costly 
noble metals (e.g., Rh, Pt) that need to be used in functioning catalysts. Moreover, the 
ability of the support material to contribute oxygen to chemical reactions can significantly 
enhance many catalytic mechanisms. Zirconium dioxide (or zirconia) is a technologically 
important catalytic support medium. This material has many other applications, as well, 
including gas sensors, solid fuel cells, and high durability coatings. 

Copper has recently been identified as an important catalytic agent for NO reduction and 
CO oxidation. 00 For instance, Cu-ion exchanged zeolites such as Cu/ZSM- alumina- 
supported CuCi and Zr02-supported coppemu all exhibit high activity for the catalytic 
promotion of NO x reduction by CO (forming N 2 and CO2). An attractive quality of the 
Cu/Zr02 system is its unusually high catalytic activity for the NO-CO reaction at very low 
temperatures (100-200°C) J0 

Because of the interest in Cu/Zr02 as a part of the next generation of automotive 
catalysts, we have performed a series of ab initio calculations to study the energetics of 
molecular adsorption on this surface. Although there have been previous theoretical studies 
of bulk ZrOolHil and a recent thorough investigation of Zr02 surfaces^!, this investigation 
is the first to study molecular adsorption onto metal films on this surface. To the best of 
our knowledge, this is also the first ab initio study of oxide-supported metal chemisorption. 
A monolayer of copper was chosen for several reasons. The monolayer geometry can be 
conveniently determined theoretically, and we believe that it can be can be reproducibly 
prepared experimentally. Adsorption onto a monolayer of copper also forms the starting 
point for systematic studies of clusters of atoms on surfaces. Furthermore, interactions of 
molecules with individual atoms in a copper monolayer makes contact with zeolite systems. 
To aid in interpretation, we also make a direct comparison between adsorption onto Cu/Zr0 2 
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and onto the bare Cu(100) surface.0 



II. METHOD 



All calculations in this study use density functional theory0'0 (DFT) within the plane- 
wave pseudopotentialil method to determine all structural parameters and energies. Opti- 
mized pseudopotentials0 for all elements were constructed to be well converged for a 50 Ry 
plane wave cut-off energy 

For structural parameters, the local density approximation (LDA)El of DFT gives highly 
accurate results for most bulk and surface systems. Typically, the optimal computed pa- 
rameters are within 1% of the experimental values. However, for adsorption energies, the 
generalized gradient approximation (GGA)Ei typical yields significantly more accurate re- 
sults. To obtain GGA results for the adsorption energies of copper and carbon monoxide 
on the substrate we take a hybrid approach. First, the structure is fully relaxed within the 
LDA (i.e., the atoms are moved until all atomic forces are less than 0.01 eV/A). Second, the 
GGA energy of the optimal configuration is obtained from the ground-state charge density 
of the LDA calculation. This hybrid approach yields significantly more accurate adsorption 
energies than LDA calculations alone; for most systems, this approach gives results very 
close to fully self-consistent GGA calculations.0 

For bulk systems, Monkhorst-Packi! special fc-point sets were used for Brillouin zone 
integrations. For the bulk cubic system, calculations using 10 irreducible Appoints gave 
convergence error of less than 1 meV per unit cell. Similar convergence criteria were used 
for the other bulk systems. For calculations involving the (111) surface unit cell, we used 
the Ramirez- BohmEU Appoint sampling method and found that convergence was obtained 
with 3 irreducible fc-points. 



III. RESULTS AND DISCUSSION 
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A. Bulk Zirconia 



The ground-state phase of zirconia — baddeleyite — has a very complex monoclinic struc- 
ture (m-ZrC^) with nine internal degrees of freedom and four formula units per primitive 
cell. Each Zr cation is 7-fold coordinated by oxygen in the monoclinic phase. The struc- 
ture can be viewed as alternating layers of Zr and O with the coordination of the O atoms 
alternating between 3 and 4 from one oxygen layer to the next. 

At about 1400 K, a first-order martensitic phase transition occurs, yielding the tetragonal 
phase (t-ZrC^). The tetragonal form of zirconia can be viewed as a simple distortion of the 
cubic fluorite structure, with alternating columns of oxygen atoms along one crystallographic 
axis shifting upward or downward by an amount d z . This structure is described by two lattice 
constants, a and c, and it has two formula units per unit cell. 

Above about 2650 K, zirconia assumes the cubic fluorite structure (c-ZrC^). In this phase 
there are only two degrees of freedom: the lattice constant, a, and an internal coordinate, 
u, reflecting the positions of the oxygen atoms along the body diagonal of the cubic cell. 
For the ideal fluorite structure, the value of this coordinate is 0.25. The cubic structure 
can be stabilized at room temperature by incorporation of a few percent of Y2O3. Addition 
of cation impurities also improves the thermochemical properties substantially, giving the 
cubic phase extremely high strength and thermal-shock resistance. 

Table 1 shows the optimized structural parameters for the monoclinic, tetragonal and 
the cubic phase of bulk zirconia (see Fig. 1). All computed values are in excellent agreement 
with experiment.il 

B. Cubic Zirconia (111) Surface 

Table 2 shows the relaxation data for two (111) surfaces of c-ZrC^: the stoichiometric 
and the fully reduced (i.e., top layer of oxygen removed) surfaces. The stoichiometric surface 
shows an outward relaxation of the top oxygen-zirconium spacing and a significant reduction 
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in the oxygen-oxygen interlayer spacing between the two outermost layers of Zr0 2 formula 
units. For the reduced surface, removal of the top layer of oxygen changes these relaxations 
dramatically. We predict a huge (24.4%) inward relaxation of the top zirconium-oxygen 
spacing and a reversal in the relaxation of the oxygen-oxygen spacing between the top two 
layers of formula units. Our calculations show a large energetic cost of 9.4 eV/atom for 
completely reducing the (111) surface. The possibility of surface reconstructions has not 
been considered in this study. A detailed study of many zirconia surfaces can be found in 
reference 16. 

C. Copper on c-Zirconia 

Before investigating the chemisorption of CO on zirconia supported copper, we first 
determine the preferred binding site for copper on (111) c-Zr02- There are three high- 
symmetry positions on this surface: the top oxygen site (01), the three-fold hollow site 
directly above the zirconium atom (Zr), and the hollow site above the subsurface oxygen 
atom in the top Zr0 2 formula unit (02). Table 3 shows the GGA binding energies for a 
monolayer of copper centered on each of these three sites. As shown in the table, the 01 
site is preferred by about 0.3 eV/atom over the other possible sites. 

It is important to compare the binding energy of copper on the c-Zr02 surface with the 
cohesive energy of bulk copper. Since the cohesive energy is nearly twice the adsorption 
energy of copper on c-Zr0 2 , the oxide-supported copper monolayer is, at best, a metastable 
configuration. Annealing this monolayer would result in the formation of copper particles 
on the oxide surface. 

D. CO on Cu on Zr0 2 (lll) with comparison to CO on Cu(100) 

With the preferred binding site of copper determined, the binding energy of carbon 
monoxide on the oxide-supported monolayer of copper can now be calculated. Table 4 
shows the GGA binding energy and the relaxed structural parameters of carbon monoxide 
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on the oxide supported copper monolayer. It is illustrative to compare the binding energy 
of carbon monoxide on zirconia-supported copper to the binding energy of a half monolayer 



distance of CO molecules is very similar for these systems. 

Our computations reveal that the binding energy of CO on zirconia-supported copper 
is nearly 0.2 eV greater than that of CO on the Cu (100) surface. This is most likely due 
to charge transfer from the copper atoms to the zirconia substrate. It has been shown 



fer from the weakly anti-bonding 5a orbital of the CO molecule to the metal conduction 
bands. This process is therefore enhanced due to the presence of the strongly electronegative 
oxide substrate. 



We have computed the bulk parameters of monoclinic, tetragonal and cubic zirconia 
which are in excellent agreement with experiment. Our surface relaxation data for the 
stoichiometric (111) c-Zr02 surface shows an outward relaxation of the top layer of oxygen. 
This may be related to the propensity of the Zr0 2 top-site oxygen atoms to participate in 
surface chemical reactions. Furthermore, we have found that the preferred binding site for 
copper on the c-Zr02 (HI) surface, by 0.3 eV/atom, is atop the surface oxygen atom. We 
have calculated the adsorption energy of CO on an oxide-supported monolayer of copper. 
Comparing our results to that of CO on Cu (100), we find that the presence of the support 
increases the binding energy by over 0.2 eV/molecule. 

ACKNOWLEDGMENTS 

This work was supported by the Laboratory for Research on the Structure of Matter as 
well as NSF grant DMR 97-02514. The authors would also like to thank the Alfred P. Sloan 



of carbon monoxide on the Cu(100) surfacet-1 As 




seen in Table 4, the nearest-neighbor 




that the Cu-CO chemisorption bond is dative, resulting from charge trans- 



IV. CONCLUSIONS 



7 



Foundation for support. Computational support was provided by the National Center for 
Supercomputing Applications and the San Diego Supercomputer Center. 



8 



REFERENCES 

1 S. Recchia, C. Dossi, A. Fusi, R. Psaro, R. Ugo and G. Moretti, Chem. Commun. 19, 
1909, (1997). 

2 D. Nachtigallova, P. Nachtigall, M. Sierka and J. Sauer, Phys. Chem. Chem. Phys. 1, 
2019, (1999). 

3 D. Biglino, H. Li, R. Erickson, A. Lund, H. Yahiro and M. Shiotani, Phys. Chem. Chem. 
Phys. 1, 2887, (1999). 

4 C. Dossi, S. Recchia, A. Pozzi, A. Fusi, V. Dalsanto and G. Moretti, Phys. Chem. Chem. 
Phys. 1, 4515, (1999). 

5 G. L. Bauerle, G. R. Service and K. Nobe, Ind. Eng. Chem. Prod. Res. Dev. 11, 54, 
(1972). 

6 Y. Okamato, T. Kubota, H. Gotoh, Y. Ohto, H. Aritani, T. Tanaka and S. Yoshida, J. 
Chem. Soc, Faraday Trans. 94, 3743, (1998). 

7 Y. Okamato, H. Gotoh, Y. Ohto, H. Aritani, T. Tanaka and S. Yoshida, J. Chem. Soc, 
Faraday Trans. 93, 3879, (1997). 

8 H. J. F. Jansen, Phys. Rev. B 43, 7267, (1991). 

9 R. Orlando, C. Pisani, C. Roetti, and E. Stefanovich, Phys. Rev. B 45, 592, (1992). 

°M. Rosenbauer and H. J. F. Jansen, Phys. Rev. B 47, 16148, (1993). 

■ 1 E. V. Stefanovich, A. L. Shluger, and C. R. A. Catlow, Phys. Rev. B 49, 11560, (1994). 

2 M. Wilson, U. Schonberger, and M. W. Finnis, Phys. Rev. B 54, 9147, (1996). 

3 B. Kralik, E. K. Chang, and S. G. Louie, Phys. Rev. B 57, 7027, (1998). 

4 J. K. Dewhurst and J. E. Lowther, Phys. Rev. B 57, 741, (1998). 

5 M. S. Kahn, M. S. Islam, and D. R. Bates, J. Mater. Chem 8, 2299, (1998). 

9 



A. Christensen and E. A. Carter, Phys. Rev. B 58, 8050, (1998). 
S. P. Lewis and A. M. Rappe, J. Chem. Phys 110, 4619, (1999). 
P. Hohenberg and W. Kohn, Phys. Rev. 136 (1964) B864. 
W. Kohn and L. J. Sham, Phys. Rev. 140 (1965) A1133. 

M. C. Payne, M. P. Teter, D. C. Allan, T. A. Arias, and J. D. Joannopoulos, Rev. Mod. 
Phys. 64 (1992) 1045. 

A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. Joannopoulos, Phys. Rev. B 41 (1990) 
1227. 

J. Perdew and A. Zunger, Phys. Rev. B 23 (1981) 5048. 

J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 18. 

M. Fuhcs, M. Bockstedte, E. Pehlke and M. Schemer, Phys. Rev. B 57, 2134, (1998) 

H. J. Monkhorst and J. D. Pack, Phys. Rev. B 16, 5188 (1976). 

R. Ramirez and M. C. Bohm, Int. J. of Quant. Chem. 30, 391, (1986). 

Experimental data taken from C. J. Howard, R. J. Hill, and B. E. Riechert, Acta Crys- 
tallogr. Sect. B 44 116 (1988) and P. Aldebert and J. -P. Traverse, J. Am. Ceram. Soc. 
68, 34 (1985) as presented in table V of E. V. Stefanovich, A. L. Shluger, and C. R. A. 
Catlow, Phys. Rev. B 49, 11560, (1994). 

E. Borguet and H.-L. Dai, J. Chem. Phys. 101, 9080, (1994). 



10 



FIGURES 
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FIG. 1. Unit cells for cubic, tetragonal monoclinic Zr02- Oxygen atoms are white, zirconium 
atoms are grey. 
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FIG. 2. Top and side view of the (111) c-Zr02 surface. Copper atoms were adsorbed to the 
sites listed (01, 02 and Zr). 
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TABLES 

TABLE I. Optimized structural parameters for cubic, tetragonal and monoclinic phases of 



zircoma. 



Parameter This Work Experiment 3 

Cubic 

a 5.043 5.090 

u 0.260 0.275 

Tetragonal 

a 5.047 5.050 

c 5.127 5.182 

d z 0.043 0.0574 



Monoclinic 



a 5.226 5.317 

b 5.082 5.150 

c 5.229 5.212 

9 99.18° 99.23° 

x Zt 0.2107 0.2083 

yzr 0.2778 0.2754 

z Zr 0.0417 0.0395 

xoi 0.3303 0.3317 

yoi 0.0771 0.0700 

z i 0.3500 0.3447 

x 2 0.4824 0.4792 

y 02 0.4477 0.4496 

z 2 0.7588 0.7569 



a Reference 27 
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TABLE II. Relaxation data for the (111) surface of cubic zirconia, stoichiometric and reduced. 
Data are listed as the percent deviation from the bulk terminated interlayer spacing. 



Spacing Stoichiometric Reduced 

A -Zr +5.89 

A Zr _ -0.03 -24.4 

A -o -3.88 +4.70 

Ao-zr +1.14 +0.70 

A Zr _ +0.62 -1.63 

A -o -0.12 +0.08 



TABLE III. Binding energies of copper at various sites on the (111) surface of cubic zirconia 



Site 

01 
Zr 
02 

Bulk Cu 



Binding Energy (eV/atom) 

1.8 
1.5 
1.3 
3.5 



TABLE IV. Structural data and binding energies of CO on Cu (100) and zirconia supported 
copper 

Parameter CO-Cu-Zr0 2 CO-Cu(lOO) 

Binding energy(eV) 0.86 0.65 

CO-CO distance (A) 3.61 3.56 

C-0 bond length (A) 1.140 1.139 

Cu-CO bond length (A) 1.801 1.852 
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